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We consider real symmetric and complex Hermitian random matrices with the additional sym¬ 
metry h xy = h]\r-y,N-x- The matrix elements are independent (up to the fourfold symmetry) and 
not necessarily identically distributed. This ensemble naturally arises as the Fourier transform of a 
Gaussian orthogonal ensemble (GOE). It also occurs as the flip matrix model - an approximation of 
the two-dimensional Anderson model at small disorder. We show that the density of states converges 
to the Wigner semicircle law despite the new symmetry type. We also prove the local version of the 
semicircle law on the optimal scale. 
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1 Introduction 

In 1955, Wigner conjectured that the eigenvalues of large random matrices describe the energy levels of large 
atoms m- Therefore, the distribution of the eigenvalues of a random matrix is an interesting and often studied 
object in random matrix theory. For an N x N random matrix with eigenvalues (A,;)(Y_ 1 , let — A -1 JT =1 

be the empirical spectral measure. The celebrated Wigner semicircle law m asserts that converges to the 

semicircle law given by the density yj {4 — £ 2 ) + /(27r) in the limit that the matrix size N goes to infinity. 

The Wigner-Dyson-Gaudin-Mehta conjecture in m asserts that the distribution of the difference between 
consecutive eigenvalues of a large random matrix only depends on the symmetry type of the matrix and not on 
the distribution of the entries. This independence of the actual distribution is called universality. The proof of 
this conjecture by Erdos, Schlein, Yau and Yin in is built upon establishing a local semicircle law in the 
first step (see [9] for a review). An alternative approach was pursued by Tao and Vu in [Ti] , 

Wigner’s semicircle law can be used to compute the number of eigenvalues contained in a fixed interval for a 
large random matrix. With the help of a local semicircle law such prediction can also be made in the case of a 
variable interval size as long as it is considerably bigger than A” 1 which is the typical distance of neighbouring 
eigenvalues. A local semicircle law is most commonly proved by establishing a convergence of the Stieltjes 
transform ■= A -1 — z ) -1 °f A 4 jv to the Stieltjes transform m of Wigner’s semicircle law. Then 

an interval size of A^ 1 corresponds to showing the convergence when rj = Imz is of this order. 

One of the most general versions of a local semicircle law is presented in [5] . They suppose that the random 
matrix H = (h xy ) x , y is complex Hermitian (or real symmetric), i.e., h xy = h yx for all x and y with real-valued 
random variables h xx for all x such that (h xy ) x < y forms an independent family of centered random variables. 
Besides assuming that the variances s xy ■= E| h xy \ 2 of a row sum up to one, i.e, 

X! Sx v ~ 1 (IT) 

v 

for all x which ensures that the eigenvalues stay of order 1, the most important requirement is the independence 
of the entries (up to the symmetry constraint). 

Many works in random matrix theory start with this independence assumption. However, some naturally 
arising random matrix models do not fulfill it. An example is the Fourier transform of a Gaussian Orthogonal 
Ensemble (GOE). For an A x A matrix H = ( h xy ) xy=1 the Fourier transform H = (h p < 2 )p,gez/jvz is defined 
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through 


Kq = Jf hx v exp ( _i 77 ( P x ~ 92/) 

x,y =1 ' 

for p,q £ 'Ll NT,. If H = (h xy ) xy= i is a real symmetric matrix then H = {h pq ) Ptqe i/NZ fulfills the relations 

hpq — h q p — h- q -p — h— p ^— q 

for p,q £ Tj/NT,. If the entries of H are, in addition, centered Gaussian distributed random variables such that 
{h xy \ x < y} are independent with E h 2 xx = 2E h xy for x ^ y then the entries of H will be independent up to this 
symmetry which we call fourfold symmetry. 

Interestingly, this symmetry also arises in random matrix approximations of the Anderson model. In [TJ, it is 
argued that the fourfold symmetry with a constant diagonal - called the flip symmetry - is a good approximation 
of the two-dimensional Anderson model in the regime of small disorder (see [2] for a review on random matrix 
models of the Anderson model). 

The first local law for Wigner matrices on the optimal scale rj « TV -1 (with logarithmic corrections) in the 
bulk has been proved by Erdos, Schlein and Yau in ;Bj. In PI, Erdos, Yau and Yin proved that mjv — m is of 
the optimal order (TV??) -1 in the bulk and they could extend this result to the edges in [12]. In the more general 
case with non-identical variances and the assumption (O), a local semicircle law on the scale rj « M 1 with 
M ■= (max Ij9 Sxy)^ 1 has been established by Erdos, Yau and Yin in PH- For this case, Erdos, Knowles, Yau 
and Yin obtained the optimal order (M??) -1 of — m in [5] even at the edge. A more detailed overview of 
the historical development of the local semicircle law can be found in section 2.1 of [3] - 

Our main result is a proof of the local semicircle law for random matrices possessing the fourfold symmetry. 
Despite the different symmetry type compared to the case in [5] the limiting distribution of the empirical spectral 
measure will still be Wigner’s semicircle law. The basic structure of the proof follows [5]. The main novelty is 
that not only the diagonal elements of the Green function have to be treated separately from the offdiagonal 
ones, but elements on the counterdiagonal need to be estimated separately via a new self-consistent equation. 

We conclude this introduction with an outline of the structure of the present article. In the following section, 
we introduce our model and some notation and state our main result. In section 3, we prove that the Fourier 
transform of a GOE satisfies the assumptions of Theorem 12.31 The remaining part is devoted to the proof of 
our main result. Section 4 contains a collection of the tools used in the proof which is given in the subsequent 
section. In the appendix, we show that the fluctuation averaging holds true for the fourfold symmetry as well. 

Acknowledgement: I am very grateful to Laszlo Erdos for drawing my attention to this question, for suggesting 
the method and for numerous helpful comments during the preparation of this article. Moreover, I thank Oskari 
Ajanki and Torben Kruger for useful discussions. 


2 Main Result 

For N £ N and x,y £ Z/1VZ, let be real or complex valued random variables (in the following we drop the 
^-dependence in our notation) such that f xx is real valued, E£ X y = 0 and E|^ xy | 2 = 1 for all x,y. Moreover, 
we assume that for every p £ N there is a constant p, p such that 


EKx»r < ^ (2-i) 

for all x,y £ L/NL and TV £ N. For fixed TV £ N, the entries are supposed to be independent up to the fourfold 
symmetry ( xy = ( yx = C-y-x = C-x-y for all x,y £ Z/NL. 

For TV £ N, let S = {s xy ) x y& x/NZ be an TV x TV-matrix of nonnegative real numbers such that s xy = s yx = 
S- Vt - X = S- X: - y for all x, y and S is stochastic, i.e. , for every x we have 

s xy = 1- (2-2) 

y 

Furthermore, we assume that the TV-dependent parameter M ■= (max X:y s xy ) _1 satisfies 

TV 5 < M < TV (2.3) 

for some S > 0. Note that the first estimate is an assumption on S whereas the second bound follows from the 
definition of M and (E3- 
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Defining h xy ■= sJyC, xy we obtain the Hermitian random matrix H^ N ' ) = {hxy) X} y^z/NZ which fulfills the 
following fourfold symmetry 

foxy — hyx — h—y^—x — h—x,—y (2.4) 

because of the definition of ( xy and the conditions on S. By definition, S describes the variances of H^ N \ 

Let p denote Wigner’s semicircle law and to its Stieltjes transform, i.e., 


P(x) '■= 2 ^\/(4 — x2 )+i 


,{z): =hL 


2 V4 — x 2 


da; 


(2.5) 


for x £ R and z £ C\M. For the real and imaginary part of z £ C, we will use the abbreviations E and r), 
respectively, i.e., z = E + ir] with E, rj £ K. 

With this definition the complex valued function m(z) is the unique solution of 


m(z) 


1 


= 0 


( 2 . 6 ) 


m(z) + z 

such that lmm(z) > 0 for 77 > 0. Denoting the resolvent or Green function of H by 

G(z) ■= (H — z)- 1 

and its entries by Gij(z) for z £ C\R we obtain for the Stieltjes transform ton of the empirical spectral measure 


m N (z) = —Tr G(z). 


We use the definitions of stochastic domination and spectral domain given in [5]. 


Definition 2.1 (Stochastic Domination). Let X = [X^ N \u)\u £ U^ n \N £ N) and Y = (F■' v )(u); u, £ 
lj( N \ N £ N) be two families of nonnegative random variables for a possibly N-dependent parameter set U^ N \ 
We say that X is stochastically dominated by Y, uniformly in u, if for all e > 0 and D > 0 there is a 
N 0 (e,D) £ N such that 


sup ) 
uec/w 


X (Ar) (u) > N £ Y (n \u) 


< N 


-D 


for all N > Nq. In this case, we use the notation X -< Y. If X is a family consisting of complex valued random 
variables and |X| -< Y then we write X £ 0^(Y). 


The definition of stochastic domination implies the following estimate which is important for our arguments 

I h X y\ ri sl'y 2 < M- 1 / 2 . (2.7) 

Definition 2.2. An N-dependent family D = (D^ w ^)jvgn of subsets of the complex plane with 

D (JV) C {z = E + ir] £ C]E £ [-10,10],M -1 < r) < 10} 
for every N £ N is called a spectral domain. 

In analogy to the matrix S, we define R = ( r xy ) = (E/i 2 y )^~®. If N is odd then R is an (N — 1) x (N — 1) 

matrix, otherwise it is an (N — 2) x (N — 2) matrix. For 77 > 0, we introduce the corresponding two control 

parameters 

r s (z) := 11(1 - m 2 (*)S)- 1 ||/oo_> / oo, r R {z) := 11(1 - m 2 {z)R)-% 0 0 ^ ( 2 . 8 ) 

and their maximum r(z) = maxjrs^),F^z)} (Note that Ts is denoted by T in [5]). 

For the definition of the spectral domain underlying our estimates, we define 

( 1 ( M “ 7 M -27 ) 1 

” E '= mi " r ; m - min \iW’ r(s)«Imm(s) J fOrall " e[iS + i’', E + il0]) (29) 

for 7 £ (0,1/2) and E £ R. Then, for 7 £ (0,1/2) the spectral domain S = S( 7 ) = (S^^jveN is defined as 

S (JV) ■■= {E+ ip-\E\ < 10, 77 ^ < 77 < 10}. (2.10) 
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Note that the spectral domain S differs from the spectral domain S in [5] due to the new definition of r(z). 
Besides this difference the following main result of this article has the same form as Theorem 5.1 in [5]. 

Theorem 2.3 (Local Semicircle Law). Let H be a random matrix with the fourfold symmetry (12.411 such that 
the conditions m and mu) are fulfilled. For 7 £ (0,1/2), we have 


|G,W-4 5 m(,)H^2^h + -T ( 2 , 11 ) 

uniformly in x,y and z £ S, as well as 

I m N (z) - m(z) | -< ( 2 . 12 ) 

uniformly in z £ S. 

The proof of our main result is based on studying self-consistent equations in the same way as the proof of 
Theorem 5.1 in [5] which uses one self-consistent equation for G xx — m. However, due to the fourfold symmetry 
it is no longer possible to directly show that the entries G x - x are small as in [5j. Therefore, we introduce a 
second, new self-consistent equation for G x - x . While deriving these self-consistent equations we will see that 
the expressions G xx — m for x £ Z/ NX and G x< - x for x ^ — x are connected among each other via E|/i xa | 2 and 
E hr xal respectively. Therefore, we introduce the matrix R in an analogous fashion as S is introduced in [5]- The 
corresponding control parameters T r and T 5 will appear in our estimates in section 15.31 Whereas the latter 
control parameter is present in [5] and denoted by T in there, the matrix R and the corresponding parameter 
Tfl are new in our work. The role of T in ;5] is filled by the maximum r(^) = max{T,s( 2 ), r^(z)}. Estimates 
on T similar to the ones in >5] are collected in Lemma 14.81 and Remark 14.91 

Remark 2.4. If the random variables h xy are complex valued with E h 2 xy = 0 for all x y then r/j(z) < CTs( 2 :) 
for z £ {E + i 77 ; E £ [—10,10], rj £ (0,10]} and therefore we can replace T by Tg in (12.91) . Thus, in this case, 
our estimates hold on the spectral domain used in Theorem 5.1 in W- 

To have a shorter notation in the following arguments, we introduce the ^-dependent stochastic control 
parameters 

k d (z) ■=rnax\G xx (z) - m(z)\, A g {z) ■= max \G xy (z)\, A-(.z) : = max.\G x ^ x (z)\, 

x x^y^—x x^—x 

A 0 (z) :=max{A s (z),A_(z)}, A(z) = max{A d (z),A 0 (jz)|. (2.13) 

Compared to [5] we added the control parameter A_ since the off-diagonal terms G X) _ x will be estimated 
differently than the generic off-diagonal terms. 


3 Fourier Transform of Random Matrices 

In this section, we give an example of a random matrix satisfying the conditions of Theorem 12.31 namely the 
Fourier transform (in the following sense) of a Gaussian orthogonal ensemble. 

Definition 3.1 (Fourier Transform). Let H = (h xy )^ y=1 be an N x N matrix. The Fourier transform H = 
(f l Pq ) p qe i/Nz the N x N matrix whose entries are given by 

h P q = h *y exp f' 1 77 ~ q y"> 

x,y =1 x 

for p,q £ "L/NTj. 

In the next Lemma we collect the basic properties of the Fourier transform of a Gaussian orthogonal ensemble 
which will imply the conditions of Theorem 12.31 

Lemma 3.2. Let H be a GOE and H its Fourier transform. Then the entries h pq and h rs are independent if 
and only if 

ip, d) £ {(A s), (s, r), (-r, -s), (-s, -r)}. 
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Moreover, H satisfies the fourfold symmetry (12.411 for all p,q £ Z/NZ. We have 


E\h pq \ 2 = N-\ Eh 2 pr =0 


(3.1) 


for all q and p r. 

Proof. To prove the if-part it suffices to show that H satisfies ( 1 P 1 ) which is a direct consequence of the fact 
that H is symmetric. 

Since h pq and h rs are jointly normally distributed and E h pq = E h rs = 0, it suffices to prove that E h pq h rs = 0 
and Eh pq h rs = 0 in order to show that these random variables are independent. The formula Eh Xiyi h X2V2 = 
JV _ 1 ( 4 ii 2 i 5 ! /i ! /2 + Sx iy2 5y lX2 ) together with 

( 2ir \ f N, m = 0, 

> exp —l — mx = < 

“ V N J ^ 0 , otherwise 

for m £ Z/NZ yields Eh pq h rs = TV -1 for (p, q ) £ {(s, r), (—r, — s)} and E h pq h rs = 0 otherwise. Thus, E h pq h rs ^ 
0 if and only if ( p , q) £ {(s, r), (—r, — s)}. In particular, E h 2 q = 0 for p ^ q. 

The relation h rs = h sr implies the first part of (13.11) and concludes the proof of the only-if part. □ 

Therefore, the Fourier transform of a Gaussian orthogonal ensemble fulfills all requirements of Theorem 12.31 
with s pq — N- 1 and ( pq ■■= iV _ 1 // 2 ft. pg . Because of the first result in the condition (12.21) is fulfilled. By 
the second part of m Remark [ 2 ~H is applicable. Thus, the local semicircle law holds true for these random 
matrices. 


4 Tools 

In this section, we collect the tools for the proof of Theorem 12.31 We start with listing some resolvent identities 
which are the basic tool for all our estimates as they encode the dependences between diagonal and off-diagonal 
entries of the resolvents. Computing the partial expectation of certain terms in expansions of the resolvent 
entries with respect to a minor will be an important step to derive the self-consistent equations. Thus, we 
introduce some notation in the second subsection. We conclude with the fluctuation averaging, an important 
mechanism to improve some bounds, and some estimates on m and T which are frequently used in our proofs. 


4.1 Minors and Resolvent Identities 

Let H = (h xy ) X y£%/ jvz be a Hermitian matrix and T C Z/NZ. 

Definition 4.1. We define the N x TV matrix and its resolvent or Green function G *- 1 '- 1 through 

:= 1 (i i T)l(j i T )hij, G^(z) := (H<n - z)- 1 
for i,j £ ZfNZ and for z £ <C\R. We denote the entries of G^ J \z) by G^\z). We set 

(T) 

£==£■ 

i t;ig T 

In both cases, we write (ai,..., a„, T) for ({ai,..., a n } U T). 

Note that is still a Hermitian TV x N matrix, in particular G^ exists. To estimate the resolvent entries 
we make essential use of the following relations. 

Lemma 4.2 (Resolvent Identities). For i,j , k (f T, the following statements hold: 

i (T ’ i) 

(t) = ha — z — TiiaG^ b ^ hbi. (4-1) 

a,b 


5 



Ifi,j ^ k then 


1 


1 


^(T) r (T) 
^(T) _ r (T,k) u ik u kj 

U ij ^ij ^(T) 


G 


kk 


g! t) 


G. 


(T.fc) 


^ik ^ki 

r m r (T,k) r (T) ■ 
U ii '-'a U kk 


If i ^ j then 


(T,i) (T ,j) 

r m _ _^(t) h r ( T >») - _r , W Y' r^’ j h 

Lr ij ~ ^ii 2—t rila(j aj ~ U jj /—/ '~' ia n ' a i' 


(4.2) 


(4.3) 


The proof of Schur’s complement formula, (ED, and the first identity in (14.21) can be found in Lemma 4.2 in 
EH and the second identity follows directly from the first one. Lemma 6.10 in [3] contains a proof of 03- 
Moreover, if rj > 0 then the spectral theorem for self-adjoint matrices yields 

£|G£ , M| a = ilmGm (2 ). (44) 

l ' 

This identity is sometimes called Ward identity. 

The functional calculus implies the following estimates on the entries of the resolvent: 

\(*ij\z)\ < ?7 _1 < M (4.5) 

for ij > 0 and all i,j& ’L/N’L. The second estimate holds if z £ D where D is a spectral domain. 


4.2 Partial Expectation 

For the partial expectation with respect to the er-algebra generated by H^ x, ~ x \ we introduce the following 
notation. 

Definition 4.3 (Partial Expectation). Let X be an integrable random variable. For x £ Z/.ZVZ we define the 
random variables E x X and F x X through 

E x X ■■= E[X\H^~% E x X ■■= X - E x X. 

The random variable E x X is called the partial expectation of X with respect to x. 

The symbols E x and F x are the analogues of Pi and Qi in [5] that were dehned by considering the minor H^\ 
Due to the fourfold symmetry column x, — x and row x, — x contain the same information, so the conditional 
expectation is taken with respect to the minor H^ x ~ x \ Notice that it may happen that x = —x, in which case 
H (x~ x ) is an (AT - i) x (N - 1) minor. 

Definition 4.4 (Independence). We say that the integrable random variable X is independent of T C Z/iVZ if 
X = E x X for all ieT. 

If Y is independent of x then F x (X)K = XY — E x (XE x Y) = F x (XT) and therefore 

EF x (X)y = EF x (XF) = E(IY) - EE x (XY) = 0. (4.6) 


4.3 Fluctuation Averaging 

Let D be a spectral domain, H satisfy the requirements of Theorem 12.31 and T a deterministic (possibly 21 - 
dependent) control parameter which satisfies 

M~ 1/2 < T < M~ c (4.7) 


for all z £ D and for some c > 0. 

The aim of the fluctuation averaging is to estimate linear combinations of the form 22k tikXk with special 
random variables Xf and a family of complex weights T = ( tik ) that satisfy 

0 < M < M-\ < 1. (4.8) 

k 
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Note that the family T may be A-dependent. Examples of such weights are given by Uk = Sifc = ¥,\hik\ 2 , 
Uk = A^ 1 or tik = Tik = Eh 2 fc . Recall that A{z) = max Xty \G xy (z) — 5 xy m(z)\ which is the basic quantity we 
want to estimate (cf. (12.131) 1. 

Theorem 4.5 (Fluctuation Averaging). Let D be a spectral domain, U/ a deterministic control parameter sat¬ 
isfying (14.71) and T = (tik) a, weight satisfying (14.81) . If A -< \f r then 


y tifcEfe—— 


-4 4r 2 , 

^ ^ kGkk 

A 4' 2 , 

^ ^ tik^kGk,—k 


k 


k^-k 


-4 >k 2 


uniformly in i and z £ D. If A -4 and T commutes with S then we have 


(4.9) 


y^tife(G fc fc - m) 
k 


~< Ts’J ' 2 


uniformly in i and z £ D. If A -4 and T commutes with R then we have 


(4.10) 


y ^ tikGk, — k 

k^-k 


A I ** 2 


(4.11) 


uniformly in i and z £ D. 

A similar result was proved in but due to the fourfold symmetry we need the third estimate in S3 
and (14.111) which were not present there. For the first estimate in (14.91) . there is the following stronger bound 
assuming that there is a stronger apriori bound on the off-diagonal terms, i.e., on A 0 (z) = max X jt y \G xy (z)\ (cf. 

(EH): 

Theorem 4.6. Let D be a spectral domain, 4/ and ifo deterministic control parameters satisfying (14.71) and 
T = {Uk) a weight satisfying (14.811 . If A -4 and A 0 -4 then 


y t ik ¥k 

k 


l 

Gkk 


a 


(4.12) 


uniformly in i and z £ D. 

The proof of Theorem 14.51 and 14.61 can be found in section 0 


4.4 Estimates on m and F 

For convenience, we list some elementary estimates from [5] which are often used in the following proofs. 
Lemma 4.7. There is a constant c > 0 such that for z £ {E + iry,E £ [—10,10], rj £ (0,10]} we have 

c<\m(z)\, |m(;j)| < 1 — cry, |m(,z)| < ry -1 , Imm(z) > cp. (4-13) 

Since T > Fg it suffices to prove the following lower bounds on T for 
Lemma 4.8. There is a constant c > 0 such that 


c<T{z), |1 — m 2 (z)\ 1 < F(^) 


for all z £ {E + iq; E £ [—10,10], ry £ (0,10]}. 

Remark 4.9. Since ||.R||.£oo_^oo < 1 the proof of Proposition A.2 in ffff yields that 


Clog A Clog A 

R ' ~ 1 — max± | | _ min{?y + E 2 ,9} 

for z £ {E + iry; —10 < E < 10, M _1 < q < 10} with 


6 = 0(z) 


if\E\<2, 

V K + Th l-^l > 2 > 


(4.14) 
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and k ■= H-E 1 ! — 2|. 


5 Proof of the Main Result 

This section contains the proof of our main result, Theorem 12.31 First, we establish the two self-consistent 
equations which will be the basis of all our estimates. In section 5.2, we bound the error terms in these self- 
consistent equations so that we can use them to prove a preliminary bound on the central quantity A (cf. 
(12.131) ) in section 5.3. Finally, we complete the proof of Theorem 12.31 in section 5.4 by iteratively improving the 
preliminary bound from the previous section. 


5.1 Self-consistent Equations 

The goal of this section is to establish the two self-consistent equations for the difference G xx — m and for the 
off-diagonal terms G x ^- x . As the matrices are indexed by elements in Z/1VZ it might happen that x = —x 
for x £ Z/1VZ, more precisely we have 0 = —0 in Z/1VZ and moreover if N is even N/2 = —N/2. Since the 
expansion of the diagonal term G xx by means of the resolvent identities is a bit different for x = — x and in this 
cases the entry G x - x is in fact a diagonal term we have to distinguish the two cases, x ^ — x and x = —x, in 
the sequel. 

Recall for the following lemma that s xa = E|/i 2 a | and r xa = E h xa . 

Lemma 5.1. For v x ■■= G xx — m we have the self-consistent equation 


y SxaVa F F x — 
z ' v x F m m 


(5.1) 


with the error term 

and the abbreviations 

G ax G 


T x = 


A*y. ■ 


hxx “ 1 “ A x Z x , x — X, 

hxx “I - A x + B x C x Y x Z x , x ^ x, 


{x-x) n (x) n {x) 


ax'-* xa 


G x 


b x .= y 


a.—x —x.a 


3xa G (x) 


C ■— (\h I 2 — « 1 h \ ' h c( x ) I u \ ' u 

'-'x — x\ °—x,xj vjr _ x — x i •*'—x,x / J ,t, xa y F r a —x ' ,l x,—x / J vjr _ x j) ll, ox 


—x,—x 

0*9 


(x,—x) 


Y x == (gL 1,_ x ) Y h xa G { Y x G^-i b h bx , 


'V (x,—x) 

2-^a,b 


Z x ■— 


h xa G^h bx 


F, 


7 s^(x,—x), 

’f'xa^ai, ,lj bx 


X = —x, 
, X ^ —x. 


(5.2) 

(5.3) 

(5.4) 


The self-consistent equation for G x - x is given by 


G Xj — x — m ^ ^ T'xaG a ,—a “1“ &X-) 
a^—a 


for where we defined £ x ■= £ x F £ x — £ x — £f with the error terms 


(5.5) 


£ x ■= - m 2 Y r xaG a -a + TR 2 Y r xaG aa + 

aG {x,—x} a——a 

(x,—x) 

£i -=g xx g^i_ x y k«Gir a) ^,-/ 


— m 


(x,—x) 


E r C 1 — C u 

1 xa^F a,—a — x llj x,—xi 


(x,—x) 


(x,—x) 


— sy(x) \ ' r C C \ ' r r<\ x ) n\ x ) 

X,—x / ; 1 xa'^ r axYJ X ,—a, ~ ^xx / J 1 xa\Ja.—x^*—x.-a' 


The self-consistent equation (15.11) has the same form as (5.9) in [5] and it is proved in a similar way by 
expanding by means of Schur’s complement formula and computing the partial expectation of a term in this 













expansion. However, we had to replace Pi by E x to derive it and the error term T x contains terms which did 
not appear in (5.8) from [5], (If x = —x then T x has the same form as in [5].) The term A x is exactly the same 
as A, in (5.8) of [5]. The term Z x is the analogue of Z t in [5] but the terms B x , C x and Y x are completely new 
and will require new estimates. 

The self-consistent equation (15.51) is new and does not have a counterpart in [5]. Due to the fourfold symmetry 
there is the necessity to introduce it since in contrast to the symmetry studied in [S] proving directly that the 
off-diagonal elements G x - x are small is not possible. 

As deriving this self-consistent equation follows the same line as the proof of (15.11) expanding and com¬ 
puting the partial expectation of a term in this expansion - it is not surprising that some error terms in EH 
have counterparts in (15.11) . Namely, £ x is the counterpart of Z x . Moreover, £% and £* are the error terms 
corresponding to A x and B x , respectively. 


Proof. We start with the proof of ED- For x = —x the derivation of m follows exactly as (5.9) in section 5.1 
of [Sj since E x and F x agree with P x and Q x respectively in this case. Similarly, for x ^ —x the self-consistent 
equation m will be obtained from Schur’s complement formula m with T = 0. In this case, its last term 
can be written in the form 


(x) (x,—x) (x,—x) 

— h h h 


E h —h n\ x ) h _i_ \ ' u r^\ x ) u \ ' u ri\ x ) u 

ll xa KJ a \ ) ,l bx — lb x,—x_ x _ x i t—x,x i / J /i, a;a u a - x' L —x,x i / j >^x,— x' jr _ rc \) ,l bx 
a,b a b 

+ E /ixaGir x) ^x + (G^ i -x )" 1 E h xa G^ x G%h bx 


a,b 


(5.6) 


by applying the resolvent identity (TOl) . Since the random variables h xa and are independent of H^ x ’ x ) 


we have E x 


n (x~ x ) h 

A*ab nbx 


= SxaG^aa ^ Sab- Thus, 


(x,— x) 

E ^ 

a,b 


7 ^y(x,—x) 7 

• lj xa Kjr a i ) ,l bx 


(x,—x) 

E q{x-x) 

^xa^aa 


G nx G 


(x, x) /^l( x ) /^l( x ) 

V^T n - rp v_T. 


\ ^ ^ '- T ax' jr xa ^ r^\ x ) \ ^ 

S X a'~ r aa ~ / J S xa ^ s —x,x'~ r -x,—x ~ / y 


a,—x ^ — x.a 


’ Xa r (x) 


where we used in the second step the resolvent identity (14.21) twice. By splitting the fourth summand on the 
right-hand side of (15.6j) according to E x + F x = 1, we get 


(x,—x) (x,—x) (x,—x) 

E h xa G { ^~ x) h bx = E E - [hxaG { a x b ~ x) h bx \ + E F - [hxaG^h^ 

a,b a,b a,b 

— ^ ^ SxaGaa A x S— XjX G_ x _ x B x “h Z x . 


Therefore, the results of (15.61) and (15.71) allow us to write (14.11) in the form 


G x 


= —z — m ■ 


T --£ 


SxaVa 


(5.7) 


which implies m using (EH). 

We fix x ^ —x. To derive (15.51) we apply the resolvent identity (14.31) twice to get 

(x,—x) 

(~l( x ) U I Si S^l( x ) \ ^ U ^y(x, — x) 7 

x,—x ^xx'J-x^x’^-a: i y-T X x'- T — Xj — X y ^ ' t xa'- 7 a 5 '^b,—X‘ 

a,b 


(5.8) 


Since E x h xa G (x ’- x) h b ,- x = G {x L a x) r xa S b ^ a splitting up the sum in the second term in (EH according to 
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E x + F x = 1 yields 


(x,—x) 


C — —C 1 h _i _ n r^\ x ) \ ' r r< I r2 _ <?3 _ c 

'~ Jr x,—x V - T x:r'- T _ x,—x x "t ^xx'-* — x,— x / j , xa y - T a,—a \ 


(5.9) 


where we used the resolvent identity (14.21) twice. We obtain (15.51) by adding and substracting to 2 J2 a r xaG a ,-a 
to the right-hand side of (Oil . □ 


5.2 Auxiliary Estimates 

The next lemma contains bounds on the resolvent entries of minors of H if there exists an apriori bound on A 
(Recall its definition in (12.131) 1. We will use a deterministic (possibly z-dependent) parameter 4/ which fulfills 


cM“5 < 'F < M~ c 


(5.10) 


for some c > 0 and all large enough N. 

Lemma 5.2. Let D be a spectral domain and ip the indicator function of a (possibly z-dependent) event. Let 
T be a deterministic control parameter satisfying (15.101) . If pA -< 4/ and T C N is a fixed finite subset then 

| -< vK -<: 4>, p\G^ ] \ -< 1, —-< 1, y\G\-m\-< pA, plmGf Ini to + A 

I Gu | 

uniformly in z £ D and in i,j for i ^ j and i,j T. 

Proof. This result follows by induction on the size of T using (14.131) and (14.21) . □ 

Using this result we will establish the first bounds on the error terms in the self-consistent equations in the 
next lemma. When applying the first part of the following lemma the indicator p will be defined precisely in 
such way that the condition pA -< M~ c holds, i.e., to ensure that ip A is small. 

Lemma 5.3. Let D be a spectral domain. 

(i) If ip is an indicator function such that ip A -< M~ c (for some c> 0) then 


p(A g + \A X \ + \B X \ + 10,1 + \Y X \ + \Z X \) -< pA~ + y ——’ (5-11) 

V>{\£l\ + \&x\ + I4 3 | + I4 4 |) <M 2 + — (5.12) 

uniformly in x and z £ D. 

(ii) For fixed r/ > 0 we have the estimates 

A_ < ?r 2 A_ + 2?r 3 A^ + e (5.13) 

with e -< uniformly in z £ £ C; Imu; = 77}, and 

A g ^M~ 1/2 + A_, (5.14) 

\A X \ + \B X \ + \C X \ + \Y X \ + \Z X \ ~< M 1 / 2 + A 0 (5.15) 

uniformly in x and in z G {w € C; Imw = 77}. 

Proof. In this proof we will occasionally split the index set of a summation into the parts {a ^ —a} and 
{a = —a} and use that the latter set contains at most two elements. 

In the following proof of the first part Lemma 15.21 will be applied several times with 4' = M~ c . Note that 
M -1 / 2 -< yj (Im to + A)/ (Mif) because of the fourth estimate in (14.131) . First, we assume x / —x. Applying 
the second estimate in (12771) and ( 1 2 . 21 ) to the definition of /4 x in ( 1 -4. 2 1 ) yields 


O) 

p\A x | S xx |G(c X | -(- ^ ' S X gip 


\G X gG a 
\G XX \ 


-< M - 1 + ip A 2 . 


(5.16) 
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Similarly, using the first estimate in Lemma 1721 we get ip\B x \ -< ipA 2 a . 
The representation 


n — \h \2sx( x ) _ x u _ u 

'-'X \ ll/ x,—x\ ^ —x ,—X °—X,X yjx —x,—X /-i l0 — X,X ,l x,— 


(x) 


G x .-x 


G-x 


G x 


c ’ 

'Gf XX 


which follows from the resolvent identity (14.31) . together with (12.71) implies 

¥> 10*1 -< M" 1 / 2 . 


(5.17) 


(5.18) 


To estimate Y x we need the following two auxiliary bounds: We have 


V? 


(x,—x) 

E j2 Mx,-x) 
,l 'xa ( ~ r a,—a 


a 


(x,—x) (x,—x) 

< 22 \hxa\ 2 ^G^’-a^ \ + E \hxa\ 2 <p\G2a ^ \ ~< pA a + M \ 
a^—a a=—a 


(5.19) 


where we used (12771) and m in last step. Now, we use the quadratic Large Deviation Bounds from [5] after 
conditioning on G( x ~ x \ By applying (C.4) in [7 with X k = ( x k and a k i = s 2k G^22 we § et 


¥> 


(x,—x) 


E h xk G { 2crh xl 


k^l 


, ( 'v° , r (x,-x)|2 , P (J V^° T r (x,-x) Im m + A 

-< 2^ s xk s xW \G k _ l I -< J]~ Y s ^ ImG kk ■< - 

k^l k 


Mr] ’ 


(5.20) 


where we used the second estimate in (1271) and (H71) in the second step. Thus, the representation 


\r _ ri{ x ) 

1 X — ^-x.-x 


' (x,—x) 


‘ (x,— x) 


22 h* a G ( a r X) h k ,-z I I E b- x ,iGi?‘- x) h, 


lb ,Lbx i 


a,k 


b,l 


(5.21) 


which follows from the resolvent identity (14.31) . yields (after separating the case k = —a) 


<p\Y x \ -< p 


(x,—x) 


E cot"” 


P 


(x,—x) 

E U x ) u 

%aO a _ fc n x k 

a^k 


, o Im m + A . o / Im m + A 

^^ + ^r <vA - + i^r 


. (5.22) 


Before estimating Z x we note that 


h xa G^~ x) h bx = E (\h X a\ 2 - s xa ) G^ + E h xa G {x b ’~ x) h bx , x ± -x, 


Zrn : - < 


” (x,—x) (x, — x) 

E F ’ 

a,b 

(x) (x) (x) 

= E (l ^“| 2 - s -) G aa +E /l “ G i“ ) ^ 

„ a,b a a^b 


(x,—x) 

E 

a^b 


We fix x ^ —x and apply (C.4) in [5] with Xi = Q xi and = s x { 2 G\ x ’ x ^ s^ 2 to get 

(x,—x) 


<P 


E k.ME X) h,- 

i¥=j 


' (x, — x) 

-< I E s xiSjx^P\G. 


y jx 

1/2 


(x,-x )|2 
ij 


-< 


Im m + A 

Mq ’ 


(5.23) 


where the last step follows in the same way as the last step in (I5.2UD . Moreover, (C.2) in [5] with X.-, = 
(|Cxi| 2 - l)(E|Cxi| 4 - 1) -1/2 and a t = (E|Cxi| 4 - 1 ) 1/2 s xi G^~ x) implies 


P 


( X,—X ) 


E {\Ki\ 2 ~s xi )G. 


( X, — X ) 


(x,—x) 


-< E s 2 xi (E\C xi \ 4 ^lW\Gt~ X) \ 2 


(5.24) 


where we used EH), the second estimate in (1271) and E2D in the last step. Therefore, absorbing M 7 2 into 
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the second summand we get 



(x,—x) 


(x,—x) 

r 

y\Z x \ < ip 

sr h r^ x ~ x) h . 

/ J "'Xl'G’ij "■JX 

+ p 

E (im 2 -»*)<#■“” 

i 

A 


/ Im m + A 


My 


(5.25) 


If x = —x then Z x can be bounded by the right-hand side in (15.111) similarly to the previous estimate and for 
A x in exactly the same way as in (15.161) . 

To estimate the generic off-diagonal entry G xy under the assumption that all of x, —x, y , — y are different, we 
use the expansion 


(x,-x,y,—y) 

G = ^ n(- x -v)n( x ~ x -v) I h _ V' h Q( x ~ x >y-v) hi 1 -I- 

^xy — ^yy ( 1l xy / J ,L xk^ki •Hy j > 

k,l 


^x.-y^-y.y , ^x,-x^-x^y_ ^ 


G 


(~x) 

-y-y 


G—x.—a 


which follows from applying ()4.3p twice and afterwards applying the first identity in ()4.2p twice. Conditioning 
on ( 3 (*.-*.s/.-!/) and applying (C.3) in [5] with X k = £ xk , Y t = Q y and a ki = s]J^G < ^ l ~ x ’ v ~ v ^s ]^ 2 yield 




(x,—x,y,—y) 

k,l 


(x,—x,y, y) 

P S xk \G 

k,l 


(x,-x,y,-y )\2 „ ; Imm + A 

ki I s iy -< Mr] > 


(5.27) 


where the last step follows exactly as in (I5.2HD . which implies 




If x = —x or y = —y then the proof of the last statement is easier. This finishes the proof of (15.111) . 

Now, we turn to the proof of (15.121) . The trivial estimate |E/i 2 y | < ¥.\h xy \ 2 = s xy < M _1 implies that the 
first two terms in <p\£ x \ are bounded by M _1 . By (12.71) its last term is bounded by M -1 / 2 . Splitting the 
summation in the third term of <p\£ x \ into a ^ —a and a = —a and using the estimate on |Eft, xy | 2 we obtain 
<p\£ x \ -4 </>AA_ + M -1 / 2 due to ( 12 . 21 ) . (14.131) . the fourth estimate in Lemma 15121 and (12.71) . Similarly to the 
bound on the third term in p\£ x \, we get p\£ x \ -4 v?A 2 and <p\£ x \ -4 </>A 2 . To estimate £ x we calculate the partial 
expectation in its definition which yields 

(x,—x) ( x,—x) 

c2 _ (~1 (~l{ x ) \ ^ (u2 _ \ r-l( x r~ x ) \ (-1 /~i( x ) \ ' L (~l{ x r~ x )-L 

°.t — ^ \ IL xa 1 xaj ^a,—a ' ^xx^—x^x / J / *xa'~ T aj _5 > L xb- 

a a^b 

Similarly to (15.241) the first term can be bounded by M -1 . Using (I5.2UI) for the second term implies 

v\£l\ < 


, 2 , / Im to + A 


My 

which completes the proof of (15.121) . 

Finally, we prove part (ii) of Lemma 15.31 In contrast to part (i), we fix y > 0. Since constants do not matter 
in the estimates with respect to the stochastic domination we will not keep track of y in such estimates. We 
start the proof of part (ii) of Lemma 15.31 with verifying (15.151) . First, we remark that applying (12.71) . (14.41) and 
(14.51) yields 


(T) 

h xa G 


(TO 

ab 


< 


1/2 , \ 1/2 
12 1 (r ir?( T ')i 2 ) 


£<&■ 


-< 


(y ‘imG^yj 


1/2 


< r] 


-l 


(5.28) 


for arbitrary finite subsets T, T' C N. The resolvent identity (14.31) and the previous bound imply 


(■ x ) 


^ | S x xG xx | + £ Sxa\G a 


( x ) 


£/ix&g; 


-4 M + A 0 , 


(5.29) 
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where we used (EZZD and (ESI) in the second step. The estimate 


I B 


(x,—x) 

x\ < E ■' 


(x,—x) 

n (x,-x) h 
U ak n k,-x 


k 


\G^, a \ ~< M - 1 / 2 


(5.30) 


is a consequence of (C.2) in [5] with X/. = (k,-x and cik = s ]2- x G^ak (14.41) . (14.511 and (12.21) . 

Applying (15.281) to the second and third term in (15.31) and (12.71) to the first term yields \C X \ -< Af -1 / 2 . 

To estimate Y x we start from (15.211) but (15.191) is estimated differently. Using the resolvent identity (14.21) 
twice we get 


(x,—x) 


E j 2 r (x-x) 

n xk U k-k 


(x,—x) (x,—x) 

A ^2 s xk\Gk,-k\ + ^2 s xk\Gkk\ 


kji—k fc= — k 

(x-x) iMx) Mx) I (x-x) 

, „ ^k-x'-’-x ,-fcl , „ \^kx^x-k\ , A , n , r - 1/2 

f T5T A ’ 


where the last step follows similarly to (15.291) and (15.301) . Combining this with the usage of (14.51) instead of 
Lemma 15.21 in (15.201) yields \Y X \ -< Af -1 / 2 + A 0 . We get \Z X \ -< A/ -1 / 2 by similar adjustments of (15.251) . This 
completes the proof of (15.151) . 

Before proving (15.131) we show 

A s <r? _ 1 A_+e (5.31) 

with some e -< Af -1 / 4 uniformly for z € {u? S C;Imw = y}. In case all of x, —x, y and — y are different it will 
be derived from the representation in (15.261) . For the fourth term in (15.261) we obtain 


\G x ,— x G — x,y\ , 

|GL«,-x| “ " 


(-x) 

/ y n -x,a'~ r a y 


<A_? 7 _1 + y~ 2 


\-x) 

< A- ( 

1/2 


1/2 


E^r )|2 


1/2 


(-x) 

^ ' (|^-— x,a\ S—x 7 a) 


(5.32) 


by applying the resolvent identity (14.31) and inserting S- Xya . In the last step, we applied (14.41) and (14.51) . Note 
that similarly to (|5.24|) we conclude that the second term is dominated by Af -1 / 4 . For the third summand in 
(15.261) we use the estimate 


sy(~ x ) ri{~ x ) 
x , — y '- 7 — y, y 


(x,—x) 


(-y,-x) 


q{~ x ) h q( x ~ x ) 

^ XX / J ,L xaGJ a _y 

a 


V h n(~y ~ x ) 

/ y IL -y,a^ay 

a 

r^i~ x ) 

^ -y-y 




where we used (C.2) in [5] as in the proof of (15.3(11) for the first factor and (15.281) for the second factor. The first 
summand in (15.261) is bounded by Af -1 / 2 due to (12.71) and (14.51) . Using (14.51) instead of Lemma l5/?l in (15.271) 
yields that the second term in (15.261) is dominated by Af -1 / 2 as well. 

We denote the sum of the absolute values of the first three summands in (15.261) and the second summand in 
(15.321) by e xy and set e — sup x , i xy . Then the above considerations show e -< Af -1 / 4 in this case. If x = —x or 
y = —y then estimating G xy is easier. Thus, (15.311) follows. 

Without inserting S- Xya in (15.321) and instead using (12.71) we see that the representation (15.261) implies (15.141) . 

To prove (15.131) we assume x ^ —x and consider the expansion 

(x,—x) (x,—x) 

C — C n( x ) \ ' r c I r< (~<( x ) \ ' T c — C h I y2 _ c3 _ c4 

^x,—x — y-*xx K - T —x,—x / j a,—a i xx'G*—x,—x / j xa^a^—a V - T xx'^ x — x ^— x ll 'x,—x i <^ x '-'x' 

a^—a a——a 


Obviously, the absolute value of the first summand on the right-hand side is not bigger than y 2 A_ and 
\£f x . < ?y _ 1 A 2 . We call the sum of the second and the third term on the right-hand side £ x and obtain 

\£ x \ -< Af -1 / 2 by (12.71) . Similarly as before, we get \£ 2 \ -< Af -1 / 2 by using (14.51) instead of Lemma l5lil An 
argument in the fashion of (15.301) yields \£ x \ < Af -1 / 2 . 
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Thus, by setting e x ■— 2 77 1 e 2 + |£ 2 | + |£ 4 | + \£ x \ and using 1)5.3111 we get 

\G x ,-x\ < V~ 2 A- + r 1 ^ + \ £ x\ + \ £ x\ + \ £ l\ < V~ 2 A- + 2^ 3 A 2 _ + e x . 

Since e x -< Af -1 / 2 uniformly in x the estimate (15.131) follows from the definition e := sup,,, e x . □ 


5.3 Preliminary Bound on A 

In this section, we establish a deterministic bound on A. The proof will make essential use of the self-consistent 
equations in Lemma 1 5. II 

Proposition 5.4. We have A -< M _7 / 3 T _1 uniformly in S. 


Once we have proved the two subsequent lemmas the proof of Proposition 15.41 follows exactly as in [5]. 
Lemma 5.5. We have the estimate 1(A < Af _7 ' /4 F _1 )A M _7 / 2 r _1 uniformly in S. 


Proof. In this proof, we will use Lemma I51T1 (il several times with ip — 1(A < M 7 / 4 T 4 ). Following the proof 
of Lemma 5.4 in [5] we get 


ipA d pT s A 2 + 


/ Im m + A 
Mr) 


since || -< ipA 2 + yj (Im m + A) /M r) by (15.111) . Moreover, because of (15.111) and the first estimate in (14.141) 
we have 


ipA g -< yTs A 2 + 


Am m + A 
Mr) 


Using m we get 


^ ^ (1 Vn R x y)Gy. — y — £x 

y¥=-y 


for all x ^ —x. Inverting (1 — m 2 R) and using (15.121) yield 


pA — 


max p\G x - x \ < Tj 7 max p\£ x \ -< pT R 

X^ — X Xy£—X 



I Im m + A 


Mr) 


(5.33) 


In total, we get 


pA -< ipT A 2 


' Im m + A 


Air/ 


as in (5.18) of [5]. Employing the definitions of S and <p as in the proof of Lemma 5.4 in [5] establishes the 
claim. □ 


When estimating the off-diagonal terms G x - x in (15.331) the control parameter T R appears naturally as the 
operator norm of (1 — m 2 !?)^ 1 in the same way as Tg (which is called T in [5]) is used in [5] to bound the 
differences G xx — m. 

Lemma 5.6. We have A -< AL -1 / 2 uniformly in z £ [—10,10] + 2i. 

Proof. We use the bounds |G^p| < 1 /77 = 1/2 from (14.51) and |m| < 1 / 77 = 1/2 from the third estimate in (14.131) . 
In particular, they imply |ri x | = |Ga; X — m\ < 1 and |m _1 | > 2. 

By (15.131) with rj = 2 we have 

A.^le^M- 1 ' 2 . 

5 

Thus, (15.141) implies A g ~< AL -1 / 2 . Hence, A a -< M -1 / 2 and therefore |T^,| -< AL -1 / 2 by (15.151) . 

Following now the reasoning of the proof of Lemma 5.5 in [3] we get A -< M 4 / 2 . □ 

Proof of Proposition The maximum of the two Lipschitz-continuous functions Tg and is a Lipschitz- 
continuous function whose Lipschitz-constant is not bigger than the maximum of the original Lipschitz-constants. 
Therefore, Proposition (B22 can be proved exactly in the same way as Proposition 5.3 in [5]. □ 
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5.4 Proof of the Main Result 

In the whole section let be a deterministic control parameter satisfying 

cM~ 1/2 < 4- < Af" 7 / 3 r" 1 . (5.34) 

The following proposition states that such deterministic bound on A can always be improved. This self¬ 
improving mechanism is also present in Proposition 5.6 of [5]. 

Proposition 5.7. Let T satisfy (15.341) and fix e £ (0,7/3). If A -<; ’I' then A -< F(f P) with 


F('P) ■= M _£ 'P + 



M E 

Mr] 


Proof. We will apply the results of Lemma T5.31 (11 with tp = 1. Using (15.111) we get 

A g + | Y x | -< A 2 + 


2 /Im to + A 2 . /Imm + 'P 


-<( T>P 2 + ... 

Mr] V Mr] 

because of the first estimate in (14.141) . The self-consistent equation (15.51) for G X; _ X implies the estimate 


(5.35) 


\G X ,_ x | < \m 2 


a^—a 


(E h 2 xa )G a ,-o 


\£ x \ -< T4- 2 


/Im m + 4/ 


Mr] 


(5.36) 


which holds uniformly in x. Here, we applied the fluctuation averaging (14.111) for G x - x with t xa = E h 2 a and 
(14.131) to the first summand, |E/i 2 y | < M _1 , Lemma l52?l and (14.1411 to the second summand and (15.121) to \£ x \ 
and employed T# < T and (|4.14D afterwards. 

Starting with these estimates the reasoning in the proof of Proposition 5.6 in [5] yields 


A -< T^ 2 + 


I Im m + A 


Mr] 


The claim follows from applying Young’s inequality and the condition *P < M 7 / 3 T 1 to the right-hand side of 
the previous estimate. □ 


In the following lemma we use the notation [u] for the mean of a vector v = ( Vi)i £ C N , i.e., 



I 


Lemma 5.8. If 4* is a deterministic control parameter such that A -< *P then we have [T] £ 0^(T 2 ). 

Proof. If x ^ —x then we obtain from Schur’s complement formula (14.11) and the definition of T x 

Y x = A x + B x — s X) _ x E x gL2,— x ~ ^xY x + F X ——• (5.37) 

XX 

The fluctuation averaging (14.91) with t,k = 1/N yields [F X G~ X ] £ 0^(\P 2 ). Obviously, we have \A X \ -< IP 2 and 
\B X | -< vp 2 by Lemma l52?l Lemma Hull Lemma l52?l and (12.71) imply |s x> _ x E x G^2,- x l ^ M~ 1 < ^ 2 due to the 
first estimate in (15.341) . 

Using (15.191) and the first two steps in (I5.2UI) with ip = 1 we obtain 


(x,—x) 

h ^Gi x f~ x) hp- x 

k,l 


-< 'L. 


(5.38) 


Thus, the representation of Y x in (15.211) and the application of Lemma \o?I\ yield \Y X \ -< 4' 2 . Hence, Lemma RTTTl 
implies |E x Yj,| -< >P 2 . For x = —x the relation (15.371) without the second to fourth term on the right hand side 
and \A X \ *P 2 hold true and |[T]| -< 4t 2 follows from (15.371) . □ 
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Proposition 15.41 Proposition 15.71 and Lemma f5.8l imply Theorem 12.81 along the same lines as Proposition 5.3, 
Proposition 5.6 and Lemma 5.7 in [Sj finish the proof of Theorem 5.1 in [5]. 


6 Proof of the Fluctuation Averaging 

In this section, we verify the fluctuation averaging, i.e. Theorem 14.51 and Theorem 14.61 To this end, we transfer 
the proof of the fluctuation averaging given in [3 to our setting. We only highlight the differences due to the 
special counterdiagonal terms G Xt - x . 

We start with two preparatory lemmas. The following result is the analogue of Lemma B.l in [5] whose proof 
works in the current situation as well. Recall that E x X = E[X|i7( x > _a; )] is the expectation conditioned on the 
minor hG,~ x ) and ¥ x X = X — E x X for an integrable random variable X (cf. definition 14.II and 14.31) . 

Lemma 6.1. Let 4* be a deterministic control parameter satisfying d' > N~ c and let X(u ) be nonnegative 
random variables such that for every pgN there exists a constant c p with E[X(m) p ] < N Cp for all large N. If 
X ( u ) -< uniformly in u then 

E x X{u) n ~< V n , ¥ x X(u) n ~< 'P n , E X{u) n -< ^ n 


uniformly in u and in x. 

This Lemma will be used throughout the following arguments. The trivial condition E[X(u) p ] < N Cp will 
always be fulfilled. The following Lemma which replaces (B.5) in [5j gives an auxiliary bound for estimating 
high moments of IXlfc Uk¥kGf k | when there are bounds on A = max Xj . y |G XJ/ — S xy m\ and A 0 = ma,x x ^ y \G xy \ 

(cf. CT). 


Lemma 6.2. Let D be a spectral domain. Suppose A -< 4/ and A 0 -< 4/ 0 for some deterministic control 
parameters 4/ and which satisfy (|5.1UI) . Then for fixed p £ N we have 



uniformly in TcN, |T| < p, x ^ T U —T and z £ D. 


Proof. If x = —x then the proof of (16.111 is exactly the same as the proof of (B.5) in [5]- For x ^ — x we start 
with 63- Since x, —x ^ T we obtain as in the proof of (15.61) by using the first resolvent identity (14.21) that 


(T,x) 


(T,tc,— x) 


(T,x,—x) 


J2 hx - G< a b ’ X) h bx =C^+ E h xa GZ’ x ’- x) h bx +(G { I’ x l*) E h xa G^lG^h bx , 


where we used the definition 


(T,x,— x) (T,x,—x) 

/"i(T) z, l _i_ \ ^ l _i_ \ ' l /~i(^i x )i , 

— ll 'x,—x'~ r —x 1 —x"'— x i x ' / j llj xa { G T a, — x ,l '— x , x ' / j ll, x,—x'G*— x b ,lj bx • 

a b 


( 6 . 2 ) 


The assumptions of Lemma 16.11 are fulfilled for each term of the expansion in (16.21) by (12.71) and the second 

estimate in 63- 

Similar to the proof of (15.181) we get \C^\ -< M -1 / 2 < \I/ 0 by (15.101) . Using the first step in (15.231) and the 
argument in (15.241) we get 


(T,x,— x) 


(T,x,—x) 


(T,x,—x) 

F x E h xa G { J b ’ x ’~ x) h bx 

< 

, r (T,x,-x) h 
/ J 'f'xa^ab ll bx 

+ 

E (IM 2 - S*a) G% X ’- X ) 

a,b 


a^b 


a 


where we used that fulfills (15.101) . The estimate 


(T,ai,— x) 

\ ^ j ^(TjX,— x) i 


k,l 


-< 


(6.3) 
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which follows from adapting (15.191) and the first step in (15.201) implies 


(T,#,— x) 


(\ \ ' L. x-y(T,x) 7 

y^—x^—xJ / J n 'xa^a,—x^ jr —x,b n ' bx 
a,b 


< ^f, -< 


(6.4) 


using a similar representation as in (15.211) and Lemma 15.21 By Lemma 16.11 these estimates imply 


(T,x) 


F x ^ K a G [ J b ’ x) h bx 


a,b 


■< Vo- 


Thus, the claim is obtained by applying Schur’s complement formula (14.11) to G xx and observing that |F x (h xx — 
~)| = \h X x\ A Tf -1 / 2 < 'F 0 as h xx is independent of H^ x ~ x l and E h xx = 0. □ 


Proof of Theorem The proof is similar to the proof of Theorem 4.7 on pages 48 to 53 in [5] so we only 
describe the changes needed to transfer this proof to its version for the fourfold symmetry. 

First, we use Lemma RTT21 instead of (B.5). Moreover, we have to change some notions introduced in the proof 
of Theorem 4.7. In the middle of page 49, an equivalence relation on the set {l,...,p} is introduced which 
has to be substituted by the following equivalence relation. Starting with k := (k\,... ,k p ) £ (Z/iVZ) p and 
r, s £ {1,... ,p} we define r ~ s if and only if k r = k s or k r = —k s . As in [5j the summation over all k is 
regrouped with respect to this equivalence relation and the notion of “lone” labels has to be understood with 
respect to this equivalence relation. We use the same notation k^ for the set of summation indices corresponding 
to lone labels. Differing from the definition in we call a resolvent entry G xy with x,y T maximally expanded 
if k l U —k l C TU {x, y}. Correspondingly, we denote by A the set of monomials in the off-diagonal entries G xy 
with T C kill —k l, x ^ y and x, y £ k\T (considering k as a subset of Z/iVZ) and the inverses of diagonal 

/nrN 

entries 1/Gxx with TckiU —k^, and x £ k\T. With these alterations the algorithm can be applied as in [5]. 
In the proof of (B.15) the assertion (*) has to be replaced by 

(*) For each s £ L there exists r = t(s) £ {1,... ,p}\{s} such that the monomial A x r 
contains a resolvent entry with lower index k s or —k s . 

To prove this claim, we suppose by contradiction that there is s £ L such that A^ r does not contain k s 
and — k s as lower index for all r £ {1,... ,p}\{s}. Without loss of generality we assume s = 1. This implies 
that each resolvent entry in A x r contains k\ and — k\ as upper index since A r ar is maximally expanded for all 
r £ {2,... ,p}. Therefore, is independent of k\ as defined in Definition 14.41 Using (14.611 and proceeding as 
in [5] concludes the proof of (*). 

Following verbatim the remaining steps in the proof of Theorem 4.7 in [5| establishes the assertion of Theorem 

mn □ 


Now, we deduce Theorem 14.51 from Theorem 14.61 


Proof of Theorem E3 The first estimate in QHD follows from Theorem 14.61 directly by setting T 0 ~ T and 
using A 0 < A T 0 . 

To verify the second estimate in (021) we use the fourth estimate in Lemma 15221 which implies 

|F x GS| = |F x (GW-m)H^. (6.5) 

Now, following the proof of Theorem lf. 61 verbatim with \I/ 0 = T and replacing the usage of Lemma RTFl by (16.51) 
yield the second estimate in (14.91) . 

Similarly, the third estimate in (14.91) is proved by following the proof of Theorem 14.61 verbatim with \I/ 0 — T 
and Lemma EH] replaced by 

|F x Gi T l,| -< A 0 -< T 

for which is a consequence of Lemma 15.21 and Lemma 16.II 

Next, we establish (14.101) . We start from Schur’s complement formula (14.11) with T = 0 and use (12.61) to get 



m 


h X kG^hi x - 


( 6 . 6 ) 
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Using Lemma 15721 with <p = 1 and the first estimate in (14.131) we get 


1 1 


G X x iti 

G X x 


GxxTYI 


TO| -< \1/. 


Thus, | h xx — (^P}k} h x kG^hi x — rnj \ -< T as well. Therefore, we can expand the inverse of the right-hand side 
of (16.61) around 1 /to which yields 


hxkG^hix -mj +g x (6.7) 

with error terms g x such that \g x \ -< T 2 uniformly in x. By (15.611 . (15.71) . (15.31) and (15.41) we have for x ^ —x the 
representation 

(x) 

^2 hxkG^i hix = S X aG aa — A x — B x — S- XiX G^_ x _ x + Z x + Y X + C x + S-x, x G^_l _x- (6.8) 

k,l a 

Taking the expectation E x of (|6.7p we want to prove that 


/ (*) 

v x ■= G xx - TO = TO 2 I -h xx + ^2, 
V k,l 


E xV x = TO. 2 ^ s xa v a + fx , (6.9) 

a 

where \f x \ -< 4/ 2 uniformly in x. From (15.71) we get that the sum of the first four summands on the right-hand 
side of (EH) is H( x ’ -measurable. Therefore, it suffices to show that all summands except the first one on 
the right-hand side of EH are bounded by T 2 uniformly in x. For A x and B x this follows directly from their 
definitions in EH). Since Z x = F^A^, for some random variable X x we get E X Z X = 0. The representation 
(15.171) for C x and Lemma l54?l yield \C X \ -< M -1 + M~ x ! 2 \ 1/ -< \1/ 2 by (15.101) . The bound (16.41) with T = 0 gives 
\Y X \ Hf 2 uniformly in x. If x = — x then the argumentation in [5] can be applied. This finishes the proof of 

EH- 

Therefore, since E^, + F^, = 1 we have 

IXa := ^ ' t ax V x = / ( tax^xVx A ^ ^ t ax W x V x = 777 / ( t ax S X yVy T F a 

xx x x,y 

— ^ ^ y SaxtxyVy T Fa — ^ ^ ^ax^x H - -^a? (6.10) 

x,y x 


where we used EH) with the notation F a ■■= Y x tax{fx + F x 7; x ) in the third step and in the fourth step that 
T and S commute. Note that |F a | -<; 'F 2 uniformly in a as \Y x taxE x v x \ = \Y x taxE x G xx \ IF 2 by the second 
estimate in EH- Introducing the vectors w — (u> a ) ae z/ivz and F := (F a )aez/ivz and writing (16.101) in matrix 
form we get 

w = m 2 S-w + F. 


Inverting the last equation yields 


w= (1 -to 2 ^)-^. 


Recalling the definition (12.81) we have 

IMIoo < FsIlFlU -< r s >F 2 

since |F a | -< T 2 uniformly in a is equivalent to HFHoo -< T 2 . This proves (14.101) . 
In order to prove (14.111) it suffices to verify 


E x G x ,-x = m 2 J2 ( Vh 2 xa )G a ,-a + fx 

a^—a 


( 6 . 11 ) 


with \f x \ -< T 2 uniformly in x. Then (14.111) follows from the same reasoning as in the proof of (14.1011 with S 
replaced by R and 

1Xx — ^ ^ txaGa^ — a- 
a^—a 
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To compute the partial expectation E x G Xi - x we use the expansion 


G x ,-x 


m 2 ^(E h 2 xa )G { :;- a x) +m 2 £ h xa G%L b x) h xb + m 2 (h 2 xa - Eh 2 xa ) G%’_T a x) 

a a^b a 

(x,—x) 

+ (to 2 - G xx G^l _ x )h x - x - m 2 h x - x + {G xx G^l_ x - to 2 ) ^ h xa G^'S b x) h xb 

a^b 

+ {G xx G^l _ x - to 2 ) ^ h 2 xa G^La' > 


( 6 . 12 ) 


which follows from the resolvent identities in a similar way as (EH). 

The first summand in (16.121) is .ff( x,-x )-measurable. Using (14.21) twice and adding the two missing terms 
we obtain the first summand on the right-hand side of (16.111) . The error terms originating from the usage of 
the resolvent identities and the added terms are obviously dominated by \F 2 . The partial expectations with 
respect to of the second and the fifth term vanish. For the remaining terms we use Lemma o First, 

|m 2 — G xx G^l _ x \ ~< \F because of the triangle inequality, Lemma [5^1 and the second estimate in (14.131) . Thus, 
using (EH and EH for the fourth term, the first step in (15.201) for the sixth term and (15.191) for the seventh 
term we get that these summands are dominated by 'F 2 . Similarly to (15.241) we see that the third summand is 
dominated by *F 2 using the Large Deviation Bound (C.2) in [5] and the first estimate in Lemma 15.21 Lemma 
16.11 establishes (16.111) which finishes the proof of Theorem 14.51 □ 
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